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Amount of night lights in an area is a proxy indicator of electricity 
consumption. This is interlinked to indicators of economic growth such as 
socio-economic activities, urban population size, physical capital, incidence 
of poverty. These night lights are generated by renewable and non renewable 
energy source. In this paper the behavior of night radiance RH data was 
minutely analyzed over a period of 28 hour; Visible Infrared Imaging 
Radiometer Suite Day/Night Band (VIIRS DNB) satellite earth observation 
data were used. These 28 hours and 8936 observations time series data is 
from 2 September 2018 to 4 September 2018. The behavior of night radiance 
RH data over 122 time intervals was analyzed using box plots. It was seen 
that the arithmetic mean of RH data is more sensitive than the arithmetic 
mean of first order difference of RH data. The first order difference of night 
radiance RH was regressed on night radiance over 110 intervals of time. The 
box plot of slope and intercept of this linear regression showed the behavior 
of these regression parameters over 110 intervals of time. It is seen that the 
data are more scattered with respect to slope than with respect to intercept. 
This implies that the rate of change in RH with respect to change in time has 
more variability that the intrinsic value of RH data at the sampled point 
of time. 
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1. INTRODUCTION 

Renewable and non renewable energies are sources to night time lights. These night lights illuminate 
bridges, streets, ship fleets, buildings etc. The data on night lights can be obtained from earth observation 
data of satellites. Researchers have found a positive correlation between night time light data and several 
socio-economic variables [1], The free and inexpensive satellite time series data on night lights can thus be 
used for in-depth statistical analysis with vast interdisciplinary applications. 

This night time light provides lighting in darkness. This night radiance is also an indicator of socio 
economic activities conducted in the night. As most of socioeconomic activities at night require light, the 
intensity of night time lights and the area covered by it correlate with socioeconomic activities and economic 
development. The brightness of indoor and outdoor night time light can be seen from the space due to 
satellite imagery. In areas lacking with data on economic accounting, night time satellite data can give us an 
estimate of economic growth. All vital inputs needed for survival and economic development of a society are 
provided by energy. Energy used in night time activity plays a very important role in this dynamics. 

Changes in population density can also be analysed using data on night time lights. Power outages in 
urban areas can also be spotted through these data. These power blackouts could be due to disasters and 
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natural calamities. Increase in night time radiance implies population growth, economic growth and 
improved living standards. 

This paper is based on data products from the Visible Infrared Imaging Radiometer Suite (VIIRS) 
night time sensors. They are also called Day/Night Band, or DNB. Elvidge, Zhinzhin, Hsu and Baugh 
discussed satellite pyrometry at night and also gave background information about these VIIRS data [2], Koel 
and Simon modelled crude birth rate and maternal mortality ratio of India using night time satellite images 

[3] , Mann, Melaas and Mallick used VIIRS/DNB to measure electricity supply reliability in Maharastra India 

[4] , Yu et al. claim that VIIRS data is a useful tool for evaluating poverty at county level in China [5], 
Sharma et al. combined Moderate Resolution Imaging Spectroradiometer (MODIS)-based multispectral data 
with the Visible Infrared Imager Radiometer Suite (VIIRS)-based nighttime light (NTL) data for robust 
extraction and mapping of urban built-up areas [6], Dou, Liu, He and Hue used VIIRS data for urban land 
extraction [7], 

In this paper a detailed statistical analysis of dynamics of change in Night Radiance Data RH is 
done. Real time radiance data is monitored from 2 September 2018, 21:16 hours to 4 September 2018, 1:56 
hours. This data is day and night band radiance satellite data. This is a detailed analysis of satellite data 
monitored for 28 hours. The results obtained here help us to identify underlying trends and patterns in the 
behaviour of this data. Energy consumption dynamics from renewable and non renewable sources used in 
night light illumination can be explained with these behaviours. The results obtained here can be used in 
analysis of long term time series day and night band radiance data. 

The main contributions of this study are as follows: 1) Statistical analysis of behaviour of day and 
night band radiance real time RH data for a period of 28 hours; 2) Statistical analysis of behaviour of first 
order difference of this RH data; 3) Regression of first order difference of RH on RH data over 110 time 
intervals and analysis of behaviour of 110 values of slopes and Intercepts obtained of this 28 hour data. 

The research questions of this study are: 1) How does day and night band radiance RH behave over 
a real time interval of 28 hours? 2) How is it related to its first order difference? 3)What are the underlying 
trends and patterns in these two variables? 


2. RESEARCH METHOD 

2.1. Data 

There are 8936 observations on time series of night radiance RH data measured by variable 
Radiative Output with unit W/m 2 . On different times ranging from 21:16 hours on 2 September 2018 to 1:56 
hours on 4 Sptember 2018, there are 122 intervals of values taken by factor time. This satellite data is a 
product of National Aeronautics and Space Administration (NASA). They provide real-time imagery of our 
night time world. These imagery products are from the Visible Infrared Imaging Radiometer Suite (VIIRS) 
night time sensors. They are also called Day/Night Band, or DNB. This paper uses VIIRS DNB radiance 
data. DNB imagery has wide range of applications. Their primary purpose is to support short term weather 
predictions and disaster response community [8], A number of sources contribute to the DNB signal, 
including city lights, lightning, fishing fleet navigation lights, gas flares, lava flows, and even auroras. When 
partial to full illumination from the moon is available, reflection of this lunar illumination off of ice, snow, 
and other highly reflective surfaces enable the study of ocean and terrestrial features. 

VIIRS is a scanning radiometer onboard the Suomi National Polar Partnership (SNPP) Satellite. The 
VIIRS collects visible and infrared imagery and radiometric measurements of land, atmosphere and oceans. It 
is sensitive to 22 wavelength bands, including a DNB with 750-m resolution. The DNB is sensitive to visible 
and near-infrared wavelengths ranging from daylight down to low levels of night time radiance. The ability 
of the DNB to detect the low levels of visible light present at night makes it well suited to studying night 
lights [4, 9], VIIRS sensors have higher resolution and detect low light better than older Defence 
Meteorological Satellite Program Operational Linescan System (DMSP-OLS) system [10, 11], Letu, 
Nakajima and Nishio estimated the co 2 emission by power plants using DMSP-OLS of VIIRS data [12], 

Rybnikova and Portnov claim that VIIRS data gives more accurate results than DMSP -OLS data 
when finding correlation between incidence of breast cancer and artificial lights at night [13]. Similarly Shi et 
al. claim that VIIRS data are better in forecasting total freight traffic for China than DMSP-OLS data [14], 
Satellite data can be used to monitor various parameters related to earth’s environment. Near real time 
monitoring of ecosystem using satellite remote sensing was proposed by Verbesselt, Zeileis and Herold [15]. 

2.2. Analysis of 28 hours near real time data 
a. Box plots 

Box plots provide visual summary of data in terms of quartiles and variance. Thus it is more 
powerful way of data representation than a tabular representation. By the use of boxplots important 
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information can be communicated and absorbed in a glance. It uses the median, the approximate quartiles, 
and the lowest and highest data points to convey the level, spread, and symmetry of a distribution of data 
values [16], Boxplots can also visually represent interrelationships between several variables. 

b. Simple Linear Regression 

The dynamics of change in response variable y is studied and predicted using single regressor 
variable x. We assume that each observation y can be described by the model [17]. 

y = Po + P±x + e 

Here f} Q is the intercept 
/?! is the slope 

Suppose that we have n pairs of observations say (y t , xf), (y 2 ,x 2 ) .(JV*?])- Then these 

data may be used to estimate /? (| and /? ] . Here y is the first order difference of night radiance data and x is the 
night radiance data. 


3. RESULTS AND ANALYSIS 

The pattern of change in RH night radiance is minutely analysed. These are 8936 observations on 
Time Series. On different times ranging from 21:16 hours on 2 September 2018 to 1:56 hours on 4 Sptember 
2018, there are 122 values taken by the factor time. The behavior of mean, median, first quartile (Ql) and 
third quartile (Q3) over 122 values of factor time is refelected in Figure 1. As seen from this figure and from 
the size of the boxplot, the variance of the average value over the 122 values of time is the higher than the 
variance of the median and Ql. The behavior of values in 65% confidence interval of the mean over 122 
values of time is reflected in Figure 2. This figure also shows that the variance of the median values is much 
lower than the variance of the mean values. The aritnmetic mean is very sensitive to extreme observations in 
contrast to median and quartiles. So the size of the boxplot for sample mean over different values of time is 
larger than Medain, Ql and Q3. 

The behavior of mean, median, first quartile (Ql) and third quartile (Q3) of first order difference of 
Rh values over 122 values of time is shown in Figure 3. The behavior of these sample statistics over arange 
of 65% confidence interval is portrayed in Figure 4. The size of the boxplot corresponding to mean and 
median in Figure 3 and Figure 4 shows that for the first order difference there are very few extereme 
observations. Hence the size of boxplot is same in Figure 3 and Figure 4. It is seen that the average values of 
first order difference of night radiance time series data are very closely spread than the average value of night 
radiance data. The behavior of spread of the night radiance data over the range of 122 values of time factor is 
refelected in Figure 5. The pattern in the spread of first order difference is exhibited in Figure 6. The size of 
the boxplot corresponding to standard deviation SD is same with respect night radiance data and first order 
difference of night radiance data. 

This is standard deviation calculated over 122 values of time. The visual display of slope and 
intercept over 110 values of time is exhibited in Figure 7. This boxplot shows that while regression first order 
difference of night radiance data on the night radiance data for different values of time, the slope has more 
variance than intercept. This implies that the intrinsic value of first order difference is concentrated around 0. 
But the rate of change in first order difference for a unit change in RH has a larger spread. The probability 
density of slope and intercept is exhibited in Figure 8 and Figure 9. The probability density of mean radiance 
over different values of time factor is exhibited in Figure 10. The probability density of mean of first order 
difference night radiance data is exhibited in Figure 11. The autocorrelation and partial autocorrelation of the 
night radiance time series data is given in Table 1. All the values obtained by correlating 8936 values are 
highly significant. 
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Night Radiance RH Time Series Data 
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Figure 1. Mean, median, first quartile and third quartile of night radiance over 122 values of time 


Night Radiance RH Time Series Data 
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Figure 2. Mean, median, first quartile and third quartile of night radiance over range of 65% 

confidence interval 


First Order Difference of Night Radiance RH Time Series Data 
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Figure 3. Mean, median, first quartile and third quartile of first order difference night radiance data over 122 

values of time 
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First Order Difference of Night Radiance RH Time Series Data 



Figure 4. Mean, median, first quartile and third quartile of first order difference night radiance data over 

range of 65% confidence interval 


Night Radiance RH Time Series Data 



Figure 5. The spread of night radiance RH data in terms of standard deviation and range 


First Order Difference of Night Radiance RH Time Series Data 



Figure 6. The spread of first order difference of night radiance RH data in terms of standard deviation and 

range 
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Figure 7. The slope and intercept of linear regressions over 122 values of time 


Probability Density of Slope 
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Figure 8. The probability density of slope over 110 Figure 9. The probability density of intercept over 
values of time 110 values of time 


Probability Density of Mean 



Probability Density of Mean of First Difference 
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N = 122 Bandwidth = 0.01927 


Figure 10. The probability density of mean radiance Figure 11. The probability density of mean of first 
over 122 values of time difference radiance over 122 values of time 
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Table 1. The Autocorrelation and Partial Auto Correlation of Night Radiance RH Data 


Lag 

i 

2 

3 

4 

5 

6 

7 

Autocorrelation 

0.234 

0.173 

0.186 

0.109 

0.087 

0.092 

0.055 

Partial Autocorrelation 

0.234 

0.125 

0.130 

0.03 

0.022 

0.037 

0.033 


4. CONCLUSION 

The results obtained in this paper are not means to an end but are complete on its own right. In this 
paper behaviour of VIIRS day and night band radiance RH is minutely explored. This is based on 28 hour 
time series data from 2 September 2018 to 4 September 2018. It is seen that the arithmetic mean of night 
radiance has more spread than the median and the quartiles. This is due to the fact that arithmetic mean is 
more sensitive to extreme observations. But the mean of first order difference is less sensitive to extreme 
observations. The size of the box plot of arithmetic means and medians for first order difference are same. 
This shows that it better to handle first order difference of RH data than the RH data itself. These values have 
been obtained over 122 values of time. First order difference of night radiance time series data is regressed 
over the time series night radiance data. The intercept of this simple linear regression conducted on 110 
values of time shows the inherent value taken by the first order difference. The slope of this simple linear 
regression shows the unit change in first order difference of RH over a unit change in the night radiance RH 
time series data. From the box plot of slope and intercept, it is seen that the spread of intercept is much less 
than that of slope. This implies that the rate of change in RH with respect to time has higher variance in 
comparison to the intrinsic value of RH for the sampled point of time. The probability density of mean of 
first order difference of RH is symmetric. The autocorrelation of 8936 time series data is highly significant. 
The partial autocorrelation is also highly significant. The probability distribution of slope and intercept gives 
an idea on the application of random effects models. These models will be explored in future research. 

But this detailed statistical analysis of such data can have vast interdisciplinary applications, 
especially for countries with limited and scarce data. The effect of night lights is a proxy to electricity 
consumption; and it permeates into diverse fields ranging from population, disaster management to 
economics. The intensity of night lights in an area can be taken as an indicator of socioeconomic activities 
conducted in that area. Such lights are generated by renewable and non renewable energy sources; and they 
illuminate areas inhabited by human population. Thus night radiance satellite data can be used in analysis of 
energy consumption dynamics of such areas. The intensity of night lights is highly correlated to the 
socioeconomic activity conducted in that area. The statistical analysis of night radiance data is related with 
the statistical analysis of energy consumption. 
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